arXiv: 1502.05138v2 [hep-th] 26 Feb 2015 


Inflation with Massive Vector Fields 


Junyu Liu“ Yi Wang 6 Siyi Zhou 6 

“University of Science and Technology of China, 

Hefei, Anhui 230026, P.R.China 

^Department of Physics, The Hong Kong University of Science and Technology, 

Clear Water Bay, Kowloon, Hong Kong, P.R.China 

E-mail: junyu@mail.ustc.edu.cn, phyw@ust.hk, szhouah@ust.hk 

Abstract. We investigate the coupling between the inflaton and massive vector fields. All 
renormalizable couplings with shift symmetry of the inflaton are considered. The massive 
vector can be decomposed into a scalar mode and a divergence-free vector mode. We show 
that the former naturally interacts with the inflaton and the latter decouples at tree level. 
The model in general predicts = 0(1), while in some regions of the parameter space 

large non-Gaussianity can arise. 
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1 Introduction 

Inflation is the most promising paradigm for the study of the early universe [1-3]. Despite 
its history and success, the detailed dynamics of inflation is still unknown. Especially, we are 
not sure if inflation is driven by one field, multiple fields, or quasi-single field. In the case of 
quasi-single field inflation, the isocurvatons have mass of order H [4, 5]. The model is later 
extended to higher mass scales [6, 7]. Similar observational signatures also arise in CFT with 
anomalous operator dimensions [8]. 

In previous studies, the model of quasi-single field inflation is investigated in the case 
where the isocurvatons are scalar fields. It is interesting to extend the study to higher spin 
fields. On the one hand, on colliders, we have found a full spectrum of spin-1/2 and spin-1 
fields but only one spin-0 field. Thus it is natural to expect the existence of those higher spin 
fields at relevant scales during inflation [9-11]. On the other hand, one important motivation 
of quasi-single field inflation, supersymmetry [12], relates particles with different spins. Thus 
a systematic study of the field spectrum helps for understanding the UV completion of quasi¬ 
single field inflation. 

In the present work, we study the coupling between the inflaton and a massive vector 
field. Although we are motivated to study the case where the vector field mass is M ~ 77, 
the calculation also holds for a general mass. The mass of the vector field would either be 
fundamental, or from symmetry breaking at higher energy scales compared to inflation. 

This paper is organized as follows: In Section 2, we study the possible interaction terms 
in this model. We list all possible interaction terms and identify a subset of terms which 
are relevant for our study. In Section 3, we introduce a decomposition method and show 
that in this case, the coupling to the vector field can be reduced to the coupling to a scalar. 
In Section 4, we calculate the power spectrum of this model by rotating the field basis. At 
leading non-trivial order, the result is reproduced by an in-in formalism calculation, which 
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we present in the appendix. In Section 5, we study the non-Gaussianity of this model. In 
Section 6 , we consider the multi-vector field case. In Section 7, we conclude and discuss some 
possible future directions. 

2 The renormalizable terms 

In this section, we introduce our setup for the model involving a scalar field as the inflaton and 
a vector field as an isocurvaton. We first summarize all the possible terms in the Lagrangian 
in Table 1. Here we only consider the terms with a shift symmetry in p, and the resulting 
equation of motion is at most second order. 


dimension 

Possible terms in the Lagrangian 

d =2 

A 2 

CO 

II 

A^dpp 

II 

A 4 ,F M ,F^,pA^d M p,d M pd^p 


Table 1. This table summarizes all the possible terms in the lagrangian when we consider inflation 
with a scalar and a massive vector field. Here we require the scalar field to have a shift symmetry. 
Also, terms which lead to higher than second order equation of motion is not considered. 


To summarize our setup, we can write down the Lagrangian of our model as 


£ = 


1 


1 


1 


--d^pd^p - -F^F^ - -M A^An + mA^d^p + fipA^d^p + k(A^A^ f - V sr (p) 

(2-1) 

where V sr (p) is the usual slow-roll potential. The signature of the metric is (—1,1,1,1). M 
is the mass of the vector field and m, (3 and k are constants. 


3 Decomposition 


In this section, we decompose the massive vector field into a transverse part and a scalar 
part. We consider the transformation 


such that 


An - An T dn(T , 


D^A^ =0 . 


(3.1) 


(3.2) 


Here the a field does not have mass dimension. We can perform a field redefinition a = Ala. 
Then the lagrangian takes the form 


C = y/^g -^d^pd^p - ^F^Ffjy - ^d^ad^a + ad^ad^p 

+upd fJ 'ad ll p + 7 (d^ad^a) 2 + (terms with A and a) — V sr (p) 


(3.3) 

(3.4) 


where a = jj, 7 = ^7 and u = jj. Here the term is invariant under this 

transformation and after the transformation it does not couple to the p field, so we do not 
consider it (for the same reason, £ IJ ' ul,X F IJU Fp\ is not relevant for our purpose). We do not 
consider either the coupling between A^ and d^a, which only contributes to higher order 
corrections in the perturbation theory. 
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4 Power spectrum and non-Gaussianities 

We start from the quadratic part of the Lagrangian 

C = s/^g -^d^ipd^ip - ^d^ad^a + ad^pd^a - V sr (tp) . (4.1) 

It is convenient to rotate the fields to remove the mixing terms. Before doing that, it is 
useful to note that the value of a is constrained. We can write down the kinetic matrix of 
Lagrangian as 

z = • (42) 

In order to avoid the ghost, we should make sure that the kinetic matrix is positive-definite, 
namely 

a 2 < 1 . (4.3) 

To solve the model, one can rotate field basis as 

P = aip + bx , (4.4) 

a = cip + dx , 

where and % are the redefined fields, and the coefficients a, b, c, d are constants. 

It is convenient to require the kinetic term of the new fields to be the same as free-held 
case, with the mixing term being eliminated. After solving the equations for a, b, c and d, we 
have the relation: 

a 2 + b 2 = - - k . (4-5) 

1 — a z 

Now we compute the two point correlation function of the fluctuation of the <p held 

( dp 2 ) = a 2 (di/j 2 ) + b 2 (dx 2 ) + ab (S'tpSx) + ab (Sxdi^) , ( 4 . 6 ) 

where dip 2 is understood as <L/?ki(^)<Vk 2 (0 h ere > which is the same for other helds. Note that 
there is no interaction between d'lj’ and dy, we have 

(d^ 2 ) = (d X 2 ) = ^(2vr) 3 <5( 3 )(k 1 + k 2 ) , (4.7) 

(d'lpdx) = (dxd'i/j) = 0 . (4.8) 

Finally 

(^ 2 > = r J 7? |) ( 2^ ,3 i (3 , ( k i+k 2 ). (4.9) 

When a is small, we can expand around a = 0 to get 

(dip 2 ) = ^3 (27r) 3 (5 (3 )(ki + k 2 )(l + a 2 ) . (4.10) 

This result agrees with a direct in-in calculation which treats the mixing of p and a as 
interactions. We present the in-in calculation in the appendix. 

To relate the inhaton perturbation to the curvature perturbation £, we have to identify 
the reheating surface. If the reheating surface depends both on the inhaton and the vector 
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directions, the isocurvature fluctuation would convert into the curvature fluctuation via multi- 
brid mechanism [13, 14]. However, we note that if we choose the rotation of this form 


ip = 

a = 


1 

Vl — a 2 
a 

y/l — a 2 


^ , 

V> + X , 


(4.11) 


a natural ending mechanism of our inflationary scenario only involves the ip field (which is 
proportional to the <p field). This is because V sr depends only on ip but not y. Thus no 
matter the end of inflation is triggered by a given field value, or a given Hubble time, the 
value of ip domains the end of inflation. 

Making use of the time delay formula £ = —HSip/ip, we get 



H 2 

8ir 2 eM 2 ' 


(4.12) 


In other words, the power spectrum is formally unchanged, if e is defined via the Hubble 
parameter 


H_^M% ndv \ 2 

lP~~2~\Vdip) ’ 


(4.13) 


where at the second half of the equation we have used the slow roll equations of motion of 
the ip field. The same fact can be found by noting that, although Stp is rescaled by a factor 
of 1 /vT — a 2 , p is also rescaled by the same factor. So we get the same result when using 
£ = —Hbip/p to calculate However, in terms of the original model, we have 





(4.14) 


As a result, with the convenient ey parameter, the power spectrum is expressed as 


pc - (i ~ 

k 8ir 2 e v M 2 


(4.15) 


Now we move on to include higher order terms in the Lagrangian: 


£ = 


~\ d^d^ip 


1<9 / V<9 A4 ct + ad^ipd^a + ^{d^ad^a) 2 - V sr (tp) , 


(4.16) 


where a = jj, 7 = Aq are coupling constants (The addition Lagrangian is £4 = kA 4 ). The 
ipd^ad/jip term gives only slow roll suppressed contributions to non-Gaussianity thus we are 
not considering it here. The reason is explained at the end of this section. 

We first consider the background behavior of this Lagrangian. The classical Euler- 
Lagrangian equations give us the classical equations of motion 


3 Hp + ip + V sr (ip) = a(3 Ha + a) , 

3 Ha + ij — 3 Hap) — aip + l2^a 2 (H& + < 7 ) = 0 . 


(4.17) 

(4.18) 











With the field rotation (4.11), ignoring the cross term from the bi-quadratic coupling of 
vector and the kinetic term of x field, we have the new Lagrangian of new fields 


£ = -^V-gd^d^ - V^aVsrm + , 


where 


6 = 


1 


\/T — a 
AM 


2 ’ 


0 = 7 « 4 tf 4 . 


(4.19) 


(4.20) 

(4.21) 


Now we need to consider correlators of ifb. The calculation of this model follows from [15]. 
We denote 

C = y/=jP(X,‘t) , 


where 


P(X, ip) = X- V sr (9ip) + 40X 2 = X + DX 2 - V sr (9i/i) 

x = , 


The energy density is 


^ 4/t /J7i\ 4 / / m\ 2 

d = 40 =au(m) W(m) 


-2 


p = 3M 2 H 2 = 2 XP x ~P = X + 3DX + V sr {9^) 


Define the sound speed as 


c 2 = 


P 


,x 


P,x + 2 XP } xx 


We define another useful parameter 


1 

U = -T — 1 . 


(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 


The leading order power spectrum of this model is given by 

< H 2 


1 - 


k 8ir 2 M?Cae ' 


(4.29) 


It is very useful to define the following two parameters when we study the non-Gaussianity 
of this model 


£ = XP.x + 2 X 2 P x ,x , 

A = X 2 P X) x + -X :> Px.x,x ■ 

Define some parameters that are more general than slow-roll as 

H XP x 




H 2 H 2 


= 0(e) , 


,= We =OU), 


(4.30) 

(4.31) 

(4.32) 

(4.33) 
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s = 


(4.34) 


He « 


= ^(e 2 ) , 


' = w = ° (e) 

The bispectrum is already calculated in [15]: 

.2 1 


(4.35) 


(CkiCk 2 Ck 3 ) — (27r) 7 <!>( 3 )(ki + k 2 + k3)(P^) 3 „ „ (A\ + A c + A a + A e + + A s ) , (4.36) 

k 1 k 2 k 3 


where 


K 


K — k\ + /c 2 + , 

(4.37) 

3 A k?k%k% 

^=2<“- 2 E> K* ’ 

(4.38) 

i>j i^j 1 

(4.39) 

A 0 = A e = Afj = 0(e) , 

(4.40) 

= 0(e 2 ) . 

(4.41) 


The variables like H, c s and e are evaluated at the time when the wave number K exits the 
horizon, namely, all the thing are measured at rfc = A - 1 + 0(e) = a + 0(e)- So the 

observable contributions are 


A = 


24 (Dxy 


u2u2h,2 
K 1^2^3 


+ 


4 DX 


1 


1 + 8 DX + 12 (DX) 2 K 3 1 + 2 DX K ^ 


E + k i k i + ^ E k 


Hi 


When we use the estimator /nl, conventionally 1 


(Ck!Ck 2 Ck 3 ) = (27r) 7 <5 (3) (ki + k 2 + k 3 )(pf) ( —)/nl ^ ^3^3 3 


at k± = fc 2 = A+ Now the measurable contribution is 




80 


(Dxy 


81 1 + 8 (DX) + 12 (DX) 


2 ’ 


_ 35 35 0X 

Tvz, - - 271 + 2DX ’ 


So 


Inl — /nl + Inl — ~ 


5 L>X(21 + llO-DX) 


(4.42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 


81 (1 + 2DX){1 + 6DX) ' 

First, we assume that D > 0, namely k > 0. Now it is very easy to show that this function is 
decreasing when DX increases from zero to positive infinity. It can be seen from the following 
Figure 1. The function has a limitation value —1|| = —0.57. So we can get the range 


275 

l'“l £ 486 =°' 57 - 


(4.47) 


1 Here the sign convention differs from [15]. We use the convention by WMAP and Planck for /nl- 
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DX 


Figure 1. The function f^L with D > 0. The blue curve is the graph of function /atl, while the 
red dashed line is its minimal — 


Now we consider the limiting case. When inflation happens, we have X r\j eH 2 Mp, so 


DX ~ DtHHll )‘ 


4k / rri \ 4 

W \m) 



(4.48) 


If a = -jj- is not too small compared to 
large. In this limit 


1, We can choose M ~ H so DX can be extremely 



(4.49) 

(4.50) 

(4.51) 


So the power spectrum is 


pC ^ 
* 8?r 2 AG< 


The three point function amplitude is 

kfknkZ 1 22 

A = 2 XX +2[ -K^ k ‘ k l 

i>j 


2K 2 

*74? 


kitf + -J]fcf) , 


(4.52) 


(4.53) 


and \/nl\ achieves its maximal 


\Inl\ 


275 

486 


= 0.57 . 


(4.54) 


It is also interesting to consider the case where D < 0, which implies k < 0. In this 
case, we can also plot f^L as a function of DX. (See Figure 2.) The ghost free condition 
requires 


H _ IP.y _ 3X + 6 DX 2 
H 2 ~ H 2 ~ X + 3DX 2 + V S r{eilj) 


(4.55) 
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Consider that the slow roll potential V sr dominates the energy density of inflation, we have 


DX >~\■ 

Similarly, the absence of tachyonic instability requires 


(4.56) 


, 2 _ Px _ 1 + 2 DX 
’ s ~ P x + 2 XP xx ~ 1 + 6 DX 


(4.57) 


Thus 


DX > 



(4.58) 


In Figure. 2, we show that this model can provide a large amplitude of /tvl when DX ~ 


l 

6 ' 



Figure 2. The /atl is blowing up where D < 0. The blue curve is the graph of function Jnl, while 
the red dashed line is its asymptote DX = — |. 

Finally, we comment on the term sj—gwtpd^ad . First, we expect \w<p\ < 0(1) because 
otherwise ghost arises in the perturbations, similar to the case of |a| > 1. In the parameter 
region of \w(p\ < 0(1), this term does not significantly modify the background slow roll 
dynamics. Still, naive dimensional analysis suggest that this term would produce large non- 
Gaussianity for natural w. However, following the formalism of [15], we find that this term 
does not modify the sound speed, nor the A/S term. Thus no large non-Gaussianity is 
generated. This fact can alternatively be observed as follows: The perturbation of this term 
can be transformed to a term which is proportional to the field equation of motion 

^TXjWipd^ad^ ~ _ ^ 3/2 y/^gw'ipd^d^ ~ - „^ 2 )3/2 , (4.59) 

where we first neglected irrelevant x terms, and then performed integration by parts. Thus 
one can redefine ^ to cancel this term. The field redefinition only leads to a slow roll 
suppressed contribution to the local non-Gaussianity thus we are not considering this term 
here. 
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5 Inflation with a tower of massive vectors 

A fundamental vector field may follow from compactification of extra dimensions, or massive 
string oscillation states. In either case, one gets not only one, bnt instead a tower of massive 
states. Following the above motivation, in this section we consider a tower of massive vectors, 
which is a direct generalization of the previous sections. 

We consider the Lagrangian 


where 


N N 


£ = E4? + E £ ?+4, 

(5.1) 

i= 1 i=l 


£-2 = 44 - l -MfAfA^ , 

(5.2) 

62 = y/^gmid^A^ , 

(5.3) 

4 = V^g (-^cpd^ - I4r(<p)^ . 

(5.4) 


This is the mass tower vector set {T/'W with corresponding field strength 

{ F$ = - d v A$ \ , mass {M i } i=12 _ N and coupling constants {ra ; } i=1 2 N . 

Decompose the vectors as 

Af = Af + d,a^, (5.5) 

satisfying 

D^Af = 0 . (5.6) 

Then we redefine the a field cr® = . We have 


/’(*) 

L -i 


V=.9 


mi 

Ah 




(5.7) 


and 


£ 


(<) = 44 = cf + 4 ° 


(5.8) 


Notice that the A® modes only couples to themselves. So we can ignore them and 
effectively write our total Lagrangian as 

N r— N r— 

£ = -'52 + ^2 ~ - ^f-gVsrip) , (5.9) 

i =1 i =1 

where a t = 4. This is similar to the single vector case. Notice that we still have the 
constraint, 

N 

E«? < 1 , (5.10) 

i =1 
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in order to avoid ghosts. 



Now rotate our N + 1 scalars as 


v = i -4 = /# » 


(5.11) 




+ x w = + x w , 


(5.12) 


where ip and {%W} . =1 2 A , are new scalars, and // = , 1 The Lagrangian becomes 


£ = - ^ —- 9 / 1 x W ^*X (t) - - V^gV sr (nip) 


(5.13) 


The interaction terms of % and ip have been canceled. 

The two point correlation function of the inflaton fluctuation can be calculated as 


(fy’k/Pka) = ^3 (2vr) 3 5 (3) ( k i + k 2 )M = ^3 (27r) 3 5 (3) (ki + k 2 )- ^ — 

2 ^T 1 - £ O? 

2=1 


(5.14) 


Again, the correction in is canceled by the correction in <p and thus the power spectrum 
for the curvature perturbation is unchanged if expressed in terms of H and e. 

For the bispectrum, at the leading order, it is similar to the N = 1 case. Namely, 
the non-Gaussianities are also at 0(e) if no nonlinear interaction is added. Similar to our 

N (i) 

previous calculation, the new Lagrangian £ A upAjp still gives trivial power spectrum 

2=1 

and bispectrum, up to slow roll suppressed corrections. 

Now add interaction 

N 




(5.15) 


Decompose similarly to the single vector case, and drop the higher order mixing, we have 
effectively 


N 

2=1 1 

After rotation, and drop the \ fields again, we have 


(5.16) 


^(*« 4 (e^) • 


(5.17) 


Similarly, we write the total effective Lagrangian 


£ = ^-gP(X^) 


(5.18) 
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where 


(5.19) 


X = --{di>f , 


and 


/ N 44' 

p=x+4 


Mf 
1 


x 2 - VM). 


Now the new D parameter 


D 


tower — 


/ N 44' 
1 / /j a t K t 

Mf 

\i=l 1 , 


N 4 
A V m i Ki 

^ Mf 
1=1 1 

N o 

Mf 

i=l * 


2 ’ 


(5.20) 


(5.21) 


Now we can conclude that the mass tower still only changes the parameter D which can be 

factored out. Now we only consider the case where all k* > 0. At the leading order, with 

N 2 

choosing suitable tower we can let ^ vA ~ 1 and cause a dramatically large DX, which will 

i =1 M * 

drive the leading non-Gaussianities |/jvl| to 0.57 for positive D towei . For negative H tower , 
large non-Gaussianity can be achieved similarly when UtowerAl is near —1/6. 


6 Conclusion and discussion 

In this work, we study the coupling between massive vectors and the inflaton field. The 
model on the one hand can be considered as an extension of quasi-single field inflation, and 
on the other hand differs in the sense that the actual degree of freedom which couples to 
the inflaton sector is massless, as long as we restrict our attention to the tree level. Thus 
equilateral non-Gaussianity, instead of quasi-local family of non-Gaussianities is generated. 

We calculate the power spectrum and non-Gaussianity of this model and find that the 
leading order non-Gaussianity is in general of order one, while having the possibility to 
become large. Then we consider the multi-vector case. We found that the power spectrum 
and non-Gaussianity take a similar form to the single vector case. 

It is interesting to consider the model beyond the leading order in perturbation theory. 
Then the divergence free part of vector start to couple to the scalar in loops and may drive 
the model towards traditional quasi-single field type. It is also interesting to study massive 
fields with other spins, especially the spin-1/2 case. We shall leave them to future work. 
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A The in-in calculation of the power spectrum 

In this part we use the in-in formalism [16] to compute the two point correlation function of 
ip up to 0(a 2 ). For a recent review of the in-in formalism, see [17, 18]. 


A.l Preparation 

We consider the Lagrangian 


72 — Up + C a + 72,; , 

(A.l) 

where 

^ ^/^gd ll '6pd IJ 5p = ^a 3 8ip8<p - ^adi8pdi8<p , 

(A.2) 

= — - y/—gd^6adnda = -a 3 8&8a — -adidadida , 

(A.3) 

72, = •J zr g l -^-d ll '8pd l ,5(j = -a 3 ^-8&8ip + a^-drfadrfp . 
1V1 lVl 1V1 

(A.4) 

We define the conjugate momentum as dir^ = and 8i: a = Jji. Then we work out the 

Hamitonian in terms of dtp, 8a and Sir^ and separate them in to the kinetic part and 

interaction part. 

Now we use the picture transformation to derive the time derivative in the interaction 
picture (we use the subscript I to denote the configurations in the interaction picture) 

jf,._ _ <977 o 

Lf>1 d8ir v 
. <977o 

Sa ' = Bin, ■ 

(A.5) 

(A.6) 

So we have 

8tt v = a 3 8ip I , 

(A.7) 

8n a = a 3 8aj . 

(A.8) 


Finally, we remove the / subscript of fields and plug the inverse mapping into the 
interaction 77,;, we get the interaction Hamiltonian in the interaction picture 77/ 


77/ = 77^ + 77/ 2 , 

(A.9) 

where 

a 3 a 2 

77/! = —8p8p , 

(A.10) 

77 1 2 = a 3 a8a8p — aadi8adi8<p . 

The mode functions can be written as 

(A.ll) 

8<p k = 8<p£ + 8ip k = u k (f)a k + u*_ k (t)a'_ k , 

(A. 12) 

8a k = hcr^ + 8a~ = v k (t)b k + v*_ , 

and the commutation relation 

(A.13) 

[a P ,Oq] = (2vr) 3 (i (3) (p - q) , 

(A. 14) 
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(A.15) 


[ 6 P^q] = (27 t)'V 3) (p - q) . 

The solutions to mode functions are 

Uk(r) = n k (r) = -JL=(1 + ikr)e~ lkT , (A.16) 

where r = — It is convenient to define x = —kr. Now the mode functions are given by 


Uk(x) = v k {x) = -^L=(i - ix)e lx . (A.17) 

Using the in-in formulism, the expectation value for 5p> 2 is 
(5tp 2 ) = ^3 (2tt) 3 (5 (3) (ki + k 2 ) + 2 Im [j^ dt i (0| 5p> 2/ H h (h) |0)] 

p~b pis pis pi ~\_ 

+ I dt\ I dti (0| 'Hi 2 (ti)d(p 2 'Hi 2 (ti) |0) - 2Re[ / dti / dt 2 (0| 5(p 2 'Hi 2 (ti)'Hi 2 (t 2 ) |0)] 
J to J to J to J to 

(A.18) 

where dip 2 is understood as 5ip kl (t)5ip k2 (t). We can write 


H 2 


where 


(6<p 2 ) = ^3 (2vr) 3 <5( 3 )(k 1 + k 2 ) + v 0 + V 1 + V 2 


Vo = 2Im[/ dt! <0| S^niMlO)] , 

Jto 

Ui= f du rdt imM n h{tl m , 

Jto Jto 

V 2 = — 2 Re[ f dti r dt 2 ( 0 | Sp> 2 H h (t 1 )'H l2 {t 2 ) | 0 )] 

Jto Jto 


(A.19) 

(A.20) 
(A.21) 
(A.22) 


A.2 Calculation of Vo 

In this section we calculate 

U 0 = 2 Im[ [ dt i (0| 5<p 2 n h (ti) |0>] . (A.23) 

Jto 

Using the momentum version of Hamiltonian, we have 

Vo = Jp lm f t dtl (°l ^k 1 (t)^k 2 (t)a 3 (U)^ p (U)<5^p(U) |0) . (A.24) 

We use the Wick contraction to deal with the integration. The integration can be 
represented as two kinds of contractions because of the rearrangement of Sip kl (t) and d(p k2 (t). 
If the integration is divergent, we can use the ie prescription and the cutoff (which means we 
can choose the lower limit of integrations to be e~ Ne instead of 0, where N e is the efolding 
number in the inflationary process) to get a finite value. One can derive that 

m 2 H 2 r+oo 

V 0 = 2 M 2 fcf ( 27 r) 3 ^ 3) ( k i + k 2 ) Im J dx 1 e~ 2iXl . (A.25) 
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After ie process we have 


So the result is 


A.3 Calculation of V\ 

In this section we calculate 


r+oo 


dx\e 


—2ix\ _ 


'0 


v '»=-ill? (2 ’' )3{<3,(k ' +k2) 


(A.26) 
(A.27) 


Ci= / dt i / dhioin^t^n^io) 

Jt-0 Jto 

Use the momentum version of Hamiltonian, we have 


(A.28) 


Vi = [ dti [ dh (0| (a 6 {ti)^Sa p (ti)Scp_ p (h) - a(h)^p 2 6a p (h)S(p_ p (ti))Sff kl (t)8cp k2 (t) 


'to Jto 


M £ 


- a(t 1 )—q 2 6<jq(t 1 )8(p_ cl (t l )) |0) . 
Split it into four terms, we have 

*i = ^11 + Via + ^13 + V 14 , 

where 


(A.29) 
(A.30) 

(A.31) 


m 


2 rt 


Vll = W J dtl J dtl ^l 0 ' ( il ) <y ^p( il ) <J ^-p( il ) <y ^k 1 W 5 ¥ , kaW a ' (*l)^q(*l)^-q(*l) |0> , 

(A.32) 

m 2 ~ 

Vl2 = ~W J dti J dti{0\a 3 (ti)5&p< y ti)5ip_ p (ti)5ip kl (t)6(p ]i2 {t)a(ti)q 2 5a cl {t 1 )dip_ (l {t 1 )\0) 

(A.33) 

Vl3 = ~W [ ^^ [ dtl ^ °(^)^ 2 ^ 0 'p(^ 1 )^-p( t ' 1 )^k 1 W%jW° 3 (tl)^q(tl)^-q(tl) 10) 

(A.34) 


TO 


2 /•< 


to -tto 
t 


Vu = jp J dtl J dtl ^ a ^ tl '> P ^q(il)^-q(*l) |0) • 

(A.35) 

For Hu, we have 


m 2 H 2 

11 = 4M 2 kf 

With the ie prescription, we have 


r+oo 

(27r) 3 (5^(ki + k2) / dxie 


2ix\ 


So 


/•+00 • 

/ dx\e 2ixi = - 
/o ^ 


Vi, = ^^s(27T) 3 <S< s >(k 1 + ka) 


16 M' 2 k\ 


(A.36) 

(A.37) 

(A.38) 
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For V\ 2 - similarly we have 
m?H 2 


V\2 = — 


2 r+oo __ /<+oo /'-I , • 

o-(27r) 3 5^ 3 ' > (ki + k2) / dxie 2 * 11 / dx’i--^ 

i Jo Jo ^1 


(1 + ixi ) 2 2ixi 


4M 2 kf 

Similarly, this integration can be evaluated using the ie prescription, 


It is very easy to show that 
So we have 


Vi 2 + V13 = 2 Re V12 = — 


Vl 3 = V,2 • 
3 m 2 H 2 


8 M 2 k\ 


(2^) 3 <5( 3) (ki +k 2 


For V\ 4 , similarly 


Vu = T^l {2!rfsl3),M+k£ 




Using the ie prescription to evaluate the integrals, we have 


Vu = 


m 2 H 2 




Summing up those contributions, we obtain Vi as 

Vi = Vn + V12 + Vis + Vu 

m 2 H 2 


M2 k 3(2 I r)V)( ki + k 2 )(I + Ie“.) 


A.4 Calculation of V 2 

In this section we calculate 


v 2 = — 2 R e [ f dh r dt 2 ( 0 | Sip 2 n h (h)n l2 (t 2 ) | 0 )] 

J to J to 

Similarly, plug in the explicit form of Hamiltonian, we have 


rti 


V 2 = —2 Re / dt\ 
Jto 
,3 


m 


(A.39) 

(A.40) 

(A.41) 
(A.42) 

(A.43) 

(A.44) 


(A.45) 


(A.46) 


m 


>t 0 


dt 2 <0| d(p kl (t)5^ k2 {t)(a {t^—d&pit^d^pih) - a(ti)—p Sa p {ti)8p_p(ti)) 


777 771 

(^{t^—Sdpii^Sip-pit 2 ) - a(t 2 )—q z 8a cl (t 2 )8ip_ cl (t 2 )) |0) 


2x~ f+~\\ ln\ (A 47) 


Split it as the last section, we have 

V 2 = V 21 + V 22 + V 23 + V 24 

where 


(A.48) 


V 21 = ~ 2 J7 2 Re / dtl J ^2 (0| ^ kl (t)^k 2 (0 a3 ( i i)^p( i i)^-p( i i) a3 ( i 2)^d q (t 2 )^_q(t2) |0) , 

(A.49) 
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^22 = 2 ^ R ej dt\ j dt 2 ( 0 \ 8 (p kl (t) 6 ip k 2 (t)a 3 (ti) 5 a p (ti) 6 ip_ p (ti)a(t 2 )q 2 Sa Ci (t 2 ) 6 <p_ cl (t‘ 2 )\ 0 ) , 

(A.50) 


m 


to J to 
2 rt rti 


V23 = 2 M? Re J dtl J dt2 (°l 8 a p {ti) 8 (p_ p {t{)a' i {t 2 ) 8 & {l (t 2 )S<p_ tl (t 2 ) |0) 


(A.51) 


R 24 = -2^R e / dtl J dt 2 (0\5if ]il (t)5ip k2 (t)a(ti)p 2 Sa p {ti)Sip_ p (ti)a(t2)q 2 8a cl (t2)8^_ cl {t2)\0) 

(A.52) 

For V 21 , using the same technologies in the last section, we have 

m 2 H 2 


Vn = 


8 M 2 kf 


(27r) 3 5( 3 )(ki +k 2 


(A.53) 


Similarly, for V 22 
m 2 H 2 


V22 = 


o (27r) 3 c)( 3 )(ki + k 2 ) Re [ dx\ [ dx 2 — e 2ia:2 . (A.54) 

1 Jo Jx! x 2 


2M 2 k\ 

It is not hard to evaluate the integration 

^22 = ^|J(2vr) 3 5( 3 )(k 1 + k 2 )(-i log 2 - \ 1E + \N e + 1 ) , 

where 7 E is the Euler 7 constant, the N e is the e-folding number. 

V 23 is very similar to V 22 - One can derive that 


(A.55) 


^23 = 


2 M 2 k\ 

Evaluating the integration as 
m 2 H 2 


'2 TT 2 r- 1-00 r+00 1 

-2-3-(27r) 3 (5 (3) (ki + k 2 ) Re / dx 1 / dx 2 (l + ~^)e~ 2lX2 . (A.56) 

d k 1 Jo Jx\ x i 


we find that 


1,23 = ^ 3 ( 2 ») 3 ' 5<3 ’( k l + k 2 )()log 2 + \~, E - \n,~\) , 


^22 + F 23 = -r0 , ?{ (27r) 3 (5 (3) (ki + k 2 ) 


2 M 2 k\ 


(A.57) 

(A.58) 


F 24 = 


For V 24 , we have 
m 2 H 2 


2A/ 2 fcf (27r)35(3)(kl + k2)Re i dXl J xl dX ^ 1+ xl^ (1+ J 2) 6 ^ ' (A ' 59) 

Evaluating the integration as 


^ = ^ (2!r)V,(kl + k 2 )( §-3 e2 '' ,) ' 


(A.60) 


As a summary, we can derive that 


V 2 = V 21 + V 22 + V 23 + V 24 


(A.61) 
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A.5 The final result 

We can see that the singular terms in V\ and V 2 are canceled. Now we can see 

<9 2 ) = |S(2’ r ) 3 < 5(3) ( k l + k 2) + V 0 + V l + V 2 

= f^ ( 2^N <3)(kl + k 2 ) ( 1 + ^ ) . (A.62) 

which agrees with the leading result from field space rotation. 
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